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ABSTRACT 


The objective of the thesis ig to minimize the total 
time required for a driver to pass through the Lincoln 
Tunnel, travel to his destination in the Central Business 
District of New York City, wait, be serviced and then exit 
his destination, This is accomplished by requiring vehicles 
to delay in New Jersey before entering the tunnel, 

A general model for total time in the system as a 

“function of the waiting time in New Jersey is presented, A 
samplins technique is developed which tests the adequacy of 
that portion of model which predicts the conditions within 
New York City. 

The control technicues are extended to apply to 


ships entering harbor through the constrains of a channel, 


viii 





CHAPTER 1 


DESCRIPTION OF THE SYSTEM 


1.1, Introduction 


The central business district (C.B.D.) of New York 
City, south from 42nd Street, may be viewed as a large and 
complex commerical and industrial center, Over the years, 
activities of a similar nature have tended to locate 
together and thus somewhat separate, although not entirely 
independent, business districts have formed, 

There are numerous entrances into and exits out of 
the C,B.D. From the point of view of the vehicles travel- 
ing into New York City, the various business districts may 
be thought of as traffic receptors or traffic sinks.  Hence- 
forth in this thesis, these will be referred to as traffic 
sinks or simply as sinks. At the present time there is no 
attempt made to limit the number of vehicles that enter the 
C.B.D, and, as a result, large traffic congestions are 
common, Consequently, a driver entering New York City can 
expect long delays before concluding his business, 

The attempt of this thesis is to reduce these delays 
by requiring vehicles to wait before entering the C.B.D, 
Ihis will allow the particular traffic sink to clear out 
enough so that the total time a driver spends in New York 


L 





City will be reduced, One entrance, the Lincoln Tunnel, 


will be studied in detail. 


1.2. Business Districts in New York City 


Although no rigorous attempt has been made to 
delineate the New York City traffic sinks, the Port of New 
York Authority has supplied the author with sufficient infor- 
mation to describe roughly the boundaries as shown in Fig. 1. 

The author attempted to determine exactly the limits 
of one of these sinks. It was found that while difficult 
and time consuming, such an undertaking was not impossible 
and could be effected if sufficient personnel and time were 
available. Fig. l is time dependent; for example, districts 
12 and 13 are ín the process of relocating further north, 
Other districts can be expected to expand, contract or 
relocate as time goes on; this indicates continuous updat- 


ins. 


1.3. The Lincoln Tunnel 

Some of the entrances to the New York City C,B,D, 
are also shown in Fig. 1. The Lincoln Tunnel is a three 
tube, six lane facility operated by the Port of New York 
Authority and represents a cost of approximately $18,000,000 
per lane per mile, (9) Each tube has 2 lanes and is about 
7,650 feet long, The south tube accommodates east bound 


traffic exclusively, the north tube handles only west bound 








Key to Fig. 1. 


District Number Description 

1 Government Complex 

2 Financial District 

3 Import/export District 
4 Fish Market 

5 Electronics (wholesale and Retail) 
6 Jewelry District 

7 Residential 

8 Textile Manufacturing 
9 Offices 

10 Light Manufacturing 

11 Greenwich Village 

p Wholesale Meat Market 
I3 Produce Market 

14 Residential 

15 Mens Clothing 

16 Garment District 

17 Residential 

18 Hospital Center 

19 Residential 

20 Insurance Center 

21 Toy District 
22 Air Line District 
23 Offices 

24 Residential 
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FIGURE 1 
The Traffic Sinks of New York City 





traffic, and the center tube is switchable, During the 
morning rush hour traffic, the center tube is used to trans- 
port vehicles into New York and during the evening rush hour 
it is changed into a west bound roadway. At other times the 
central tube handles traffic flows in both directions, 
Although rather complete information on the origins 
and destinations of vehicles which use the Lincoln Tunnel 
has been collected(25) and data on the behavior of vehicles 
within the tunnel has been obtained, (26) no success has been 
realized in attempting to optimize the use of the center 
tube, (11) 
Since the adequacy of the control techniques, to 
be presented in later chapters, depends upon the traffic 
engineer's ability to control the traffic within the 
Lincoln Tunnel, some of the theories of traffic flow will 


be examined in the next section, 


1,4 General Theories of Traffic Flow 
Many authors have studied traffic flow in detail and 


only a brief review will be presented here, One relationship 
basic to the study of traffic flow is response at time (t+T) 
equals sensitivity multiplied by stimulus at time (r), (17) 
Response is usually taken as acceleration since that is the 
Variable under the driver's control. The stimulus is gener- 


ally represented as the difference in the velocities of two 
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2 
vehicles. Thus, — (x_(t+T))= sensitivity x Lëschte LEI minus 
at 


SX (©) where X represents distance traveled, For 


mathematical simplicity, sensitivity was originally taken as 
a constant. Later, however, more sophisticated models of 
sensitivity have been developed(i2) 

A fundamental diagram of traffic flow has been 


developed(l^) and is shown in Fig. 2. It is stated that 
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FIGURE 2. 
Fundamental Diagram of Traffic Flow 


such a diagram is valid only for a particular section of 
roadway at a particular time, Authors have placed different 
interpretations on this fundamental diagram, Martin(i?) 
feels that it is valid only in equilibrium or steady state 
conditions and as such represents the desires of the drivers. 


Duckstein believes that Fig. 2 is the center oí a band of 





6 
possible flowedensity combinations. Many experimenters feel 
that traffic density (k) is the best indication of the 


actual traffic conditions i?,21). 


Average velocity is not a 
good indication because small speed changes can represent 
large changes in flow, (24) Greenberg has compared traffic 
flow to that of a continuous fluid, (13) 

Shock waves appear to be a major cause of low flow 
as, (297 It has been shown that a small increase in peak 
flow rates results in dramatic reductions in traffic congese 
tion, (10) 

Authors feel that density can be used to predict 
the formation of a shock wave(^) and that high output speeds 
indicate the absence of a shock wave. Officials of the 
Tunnels and Bridges Department of the Port of New York 
Authority have demonstrated that there exists a critical 
output velocity and that traffic densities can be used to 
foresee subecritical velocities, (10) Finally, it has been 
shown that inherent instability exists in non-homogeneous 
traffic flows, 2) The velocity at which traffic actually 
travels is a compromise between the desired velocities of 
various components of a nonehomogeneous flow, Thus the 
components which are required to travel below their optimum 


velocity can react to a sudden decrease in the speed of 


leading vehicles in such a way as to create a shock wave, 





The Port of New York Authority has been and is 
presently attempting to induce small, frequent gaps in the 
traffic flow. These small spaces do not decay as rapidly 
as larger, less frequently spaced gaps. Only a small degree 
of success has been realized in implementing these goals, 
Ihe south tube of the Lincoln Tunnel is equipped with a 
closed circuit television monitoring system, but as of this 
writing, the tunnel officials only react to but do not pre- 
vent the formation of shock waves. 

The traffic sinks in New York City and the Lincoln 
Tunnel have been discussed briefly. This is necessary for 
an understanding of the objective of the thesis which is 


presented in the next section, 


1.5. Objectives of the Thesis 
Once a vehicle departs the tunnel, his problems have 


just begun, A driver will find himself fighting traffic 

all the way to his particular traffic sink and then fighting 
for a position in the sink, It is felt that by requiring 

a vehicle to wait in New Jersey before entering the tunnel 
it may be possible to allow his traffic sink to depopulate 
sufficiently so that the total time required for him to be 
serviced in New York, is actually reduced, The basic idea is 
to remove this driver from the congested area long enough 


for some vehicles in the area to depart, 





Since it is almost an impossible task to determine 
the destinations of cars and small vehicles, it was decided 
to apply this principle to large trucks and tractor trailers 
which supply and remove goods from the C,B,D, of New York 
City. 

The general model is presented in Chapter 2, Two 
basic problem areas arise: (1) sampling the sinks to 
estimate various parameters and (2) determining if the 
tunnel supplies a large enough percentage of the sink 
traffic to make effective control of the sink conditions 
possible. These problems are treated separately in Chapters 
3 and 4, 





CHAPTER 2 
MODELING THE SYSTEM 


2.1. General Model 

As a vehicle approaches the Lincoln Tunnel In New 
Jersey, it enters a large and complex traffic system, This 
system consists of the waiting areas in New Jersey, the 
tunnel, the network of roads in New York leading to the 
destination, and the traffic sink itself, Upon being sere 
viced and departing from the traffic sink, a vehicle is cone 
sidered to have exited the system, The Lincoln Tunnel may 
be thought of as an entrance way leading to 24 separate but 
not wholly independent subesystems or traffic sinks, As a 
vehicle leaves its sink it may travel into another sink 
and thereby become part of the uncontrollable input into 
another sub=systen 

Total time in the system (T) consists of six compon- 


6 
Y t. = T where: 


ents (ti) such that ifa 


1 


ES = waiting time in New Jersey, 

t, = travel time in the tunnel, 

t, = travel time from tunnel exlt to desired traffic 
sink, 


t, = waiting time in sink, 
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ts = service time in sink, 

tg = time required to exit sink after being serviced. 

The time required for a vehicle to return to its 
original destination is assumed constant and is not con- 
sidered in the model, All times are measured in minutes. 
This model is valid for any time period during the day, howe 
ever, here it will be applied to the morning rush hour 
traffic between the times t and t 


start end’ 
All times are assumed to be between these Limits 


See Fig. b. 


2.2. Components of the Model 
The individual components of T will now be investi- 


gated separately. 
2.2.1. Waiting Time in New Jersey 
Waiting time in New Jersey (t) is the independent 


variable. T will be minimized by controlling t An upper 


1° 
bound of 25 minutes has been chosen for this variable. 

2.2.2. Travel Time in the Tunnel 

It ege, Sait originally that in order to predict E 
satisfactorily it would be necessary for the model to predict 
the status of each vehicle in the tunnel at every second, 
Ihis involved estimating distance traveled, instantaneous 
velocity and acceleration and headway for every vehicle, It 
soon became obvious that this approach would be of little 


practical value because of the large amount of computer 


memory space and computing time required to predict tp. 





11 
To avoid these problems, it was decided to utilize 
the flowdensity data available for the ER A quad- 
ratic curve was fitted to this data by the method of least 
squares, and the solution of the normal equations was: 


02a + bk + ck? where 


D e flow measured in vehicles/minute, 
k = concentration measured in vehicles/mile, 
a = 0,44, b - 0,61, and c 2-0,0043, 


The predicting curve is shown in Fig. 3. 


“low (Vehicles/minute) 





O 60 70 o0 90 100 
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FIGURE 3 
Flow (8) vs Concentration (k) 
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It is also possible to observe flow as a function of 
time of day. This data was obtained and the flow between 


the times of t and t was approximated by a quadratic 


start end 


function, This curve is shown in Fig. 4, The quadratic 
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FIGURE 4. 
Flow (8) vs Time of Day (t) 


equation is 9 = d t? 4 pt + q with d = -9,2, p = 127,8 and 
q = -241, 

If it is assumed that no shock waves will appear in 
the traffic flow, then 90 remains fairly constant for any 
given time of transit through the tunnel. For any time of 
the day, 


ck? +bkK tae 6 = 0, 


-b4/b^-4c(a-0) 


k = — RÀ 
2c k 





rG 


— 2.6144. 612-4(- 0043) (, 64-9) 
2(-.0043) 





For t < t xb ; 9 will always be greater than 
start ~ ~ end 
.%4 but the expression under the radical will always be 
positive, This can be seen by substituting as for €, 
Since the model predicts traffic concentration for rush 
hour traffic, it seems reasonable to assume that the actual 


density is the largerof the two densities corresponding to 


a specific value of 6, Thus, 


k = 61 + 6 o4 = 043 en D 


2(.0043) i 


The average velocity through the tunnel (v) is the flow 
divided by the density. Thus, velocity in miles/minute 
may be written as; 


v = 9/k, 





v s 20c/(-b4^ b^ 4C(a-0)). (1) 
As expected, "o, l simply says that on the heavy concentra- 


tion half of Fig. 4, velocity increases then decreases with 


9, 
ae 
It was show that 6 = dt” + pt + q for t tart SËLL 
toand This expression is substituted into q. 1 to obtain 


velocity as a function of time: 


v = 2c(dt24pt*a)/(-bw/ b? -4c(aedtZepta-q)). 


Since this velocity is considered to be the average velocity 


through the tunnel, the expected time required for a passage 
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through the tunnel is found by dividing the tunnel length 
(L miles) by this average velocity. Thus, after inserting 


the values of the parameters, 


- LC Le! 617-5(2,0043) (, 4422, 2657-331:971) (2) 


(-2.,36t 433t-97)(2)(.,.0043) 


Since t, the time at which a vehicle enters the tunnel, is 
the time of arrival at the New Jersey entrance plus the 
waiting time in New Jersey, it is seen that t, depends on 


the control variable t: As t varies between E and 


start 


tenda» ? increases and then decreases (Fig, 4). This results 


in the unlikely conclusions that t, reaches its minimum 


2 
value when t = t*, or when 0 is at its maximum value, This 
is shown graphically in Fie. 5 in which t, Ls plotted as 
a function of t as t varies between Cr and Fond 

This erroneous result was cuased by the two follow- 
ine unrealistic assumptions: 

1, The actual concentration is always the larger 
of the two concentrations corresponding to a specific 0, 

2, No shock waves will be induced in the flow. 
Although it is true that during peak traffic the flow is on 
the heavy concentration side of Fig. 2, it is not reasonable 


to asume that the flow started on that side, As the morning 


rush hour approaches, flow increasesalong the left hand side 
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Time in Tunnel (Minutes) 





Time of Day (Hours) 


FIGURE 5. 
Time in Tunnel (t3) vs Time of Day (t) 


of Fig. 3. At some point, the flow shifts over to the right 
hand side of the curve; this can happen in one of two 
ways; 
1, Flow can reach its maximum from the left and 
then decrease down the right side of Fig. 3, 
2. As flow is approaching its maximum value, a 
shock wave can shift the flow to the heavy concentration 


side of the graph. Thus the true expression for t, as a 


function of t is: 
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Po AE Er 
t. = Li, 61-7, 61°-6(= 0043) (4649, 28°.127,86+421)), ¿e 
S (-9,2t*+127,8t-421)(2)(.0043) start=" 1 
or 
T ws en 
E Lee O -127,8t4421)) — 


E (-9,2t^4127,.8t-421) (2) (.0043) 


where + represents the time at which maximum flow is 
achieved if no shock wave occurs or the time of a shock wave, 

The solid lane in Fig, 6 shows how t, varies with t 
if no shcok wave is induced; the dotted line traces t. if a 
shock wave is encountered at t. In either case, it can be 
concluded the t, is a nonedecreasing function of t. 

If Fig. 4 were used as a guide, one would decide 
that the flow will never be on the light concentration side 
of the flowedensity curve until evening rush hour. This 
obviously is not necessarily true, An acceleration shock 
wave can occur as 9 is decreasing along the right side of 
Pig. 3. If this happens, t, will begin to decrease, Fig. 


7 shows this penomena graphically, 


2.2.3 Travel Time to the Traffic Sink 

ihe author had hoped to study t, in detail, but 
because of time limitations, this was impossible, Therefore 
t} will be taken as a constant for each traffic sink, Refere 
ences (2) and (23) among others, discuss some of the prob- 


lems that would be encountered in predicting t4. 
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FIGURE €. 
Time in Tunnel (t,) vs Time of Day (t) 
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2.2.4 Time in the Traffic Sink 

It is almost impossible to know exactly what will 
happen to a vehicle once it enters a traffic sink. Some 
vehicles will find a legitimate servicing station immediatee 
ly, some will create their own servicing stations by double 
parking or by parking in a restricted zone, and some 
vehicles will be forced to wait until a servicing station 
becomes available, It is also impossible to determine 
exactly what will happen to a vehicle once it locates and 
enters a servicing station, Some vehicles will find men 
available to unload or load cargo while other vehicles will 
be required to wait for servicers, As the congestion mounts, 
servicers will begin to interfere with each other and will 
increase service time, It becomes clear that the standard 
queueine models are not adequate to describe the situation, 

It is hypothesized that the total time spent in a 
eraffic sink (t, + ts + t.) increases with the number of 
vehicles in the sink at the tiae of arrival (N). It is 
further hypothesized that this increase is not linear but 
that the rate of change of the total time in a sink with 
respect to the number in the sink also increases with N, 

Letting w = t, + tz + to, tho following model will 
be hypothesized: E(w) = "T where a and B are parameters 


which must be estimated for each traffic sink, The possi- 


bility of using other models for w is discussed in Chapter 3, 
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For the moment, it is assumed that t, = tj = O so 
that vehicles enter their sinks directly and instantaneously 
from the New Jersey waitins areas. Total time in the sink 
will now be studied as a function of Zu 

If it is assumed that all entries into a sink occur 
instantaneously at the beginning of each minute, then a 
consequence of the model is that all the vehicles that enter 
a particular traffic sink during the same minute have the 
same expected tine in the sink Z(w). Trucks which enter a 
traffic sink during the same minute are said to comprise a 


block of vehicles. At time t, a given sink may be represent- 


ed as in rig. 8. Here N. represents the number of vehicles 


j 





FIGURE 8 
Conditions within a Traffic 
Sink, 


expected to depart at time t The expected number of 


3: 
entries per minute from all sources other than the tunnel 
is denoted by i. It is assumed that A remains constant 
for torartstsStang and that the first in first out rule is 
k 
valid for the traffic sinks, Letting N = 2 N, the 
=a 


following relationships can be devloped: 
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"T ` 
aep (NHA +1) 


pee 
FA 
rt 
l 


= 0, E(w) = 


N+2A + 
l 1, E(w) = aef t A +1) 


rt 
l) 


AB (N+(t +1)X+1) 


and in general E(w) =a for t, « t, - t. 


L 


B(N (t]+1)A+1) reaches its minimum 


The sum t, + E(w) = t,+ae 
value when t) = 0. 

For t. * € < t, <,€, = t The general expressions for 
expected time in the sink is E(w) = gee (NaN +(ty+LA +1) The 
sum Of t, and w is minimized when E BES = t. Similarly 
for L tam ty Es - t the total time in the system is 
goe (Ne Nau nb *Ni+(t1,+1)1 +1) and is minimized when t, - 

t, - t. 

It is also seen that the value of t, which minimizes 
I is such that a vehicle will arrive at the traffic sink at 
the expected time of departure of a block of vehicles and 
never between the expected times of departues of two blocks. 

Unfortunately t, #0 and t4 £ 0 and it is not poss- 
ible to control t, so that a vehlcle will enter atrafflc slnk 
just as a block of vehicles is departing. Considerins that 
vehicles must go through the tunnel before entering their 
desired sinks, it is clear that the present expression for 
w is not adequate, When a vehicle arrives at the tunnel, the 
control personnel must predict the conditions at his desired 
traffic sink for the coming 25 minutes. Then a decision can 


be made as to how long to delay the vehicle in New Jersey. 





“nis prediction must consider the vehicles bound for this 
specific sink which have already entered the tunnel but have 
not yet arrived at the sink, 

Letii be the number of these vehicles which have 
passed the control personnel but have not yet arrived at 
the sink, Then the expression for E(w) must include a term 
fM o EE)? which predicts the arrival of these vehicles. 
Let SG, t, ti) predict the dissipation of vehicles from the 
sink. The expression for E(w) is now 

E(w) = qexp(8B (N 8C, t,t] )+E(M ,t,t, )*(t, 1). *1)). 


2.3. Minimizing Total Time_in the System 


unl. Minn an to "Mr 


Combining all parts of the model yields: T(t,t,)=t,+ 


L(.61-r/,612-4(-.0043) (.6449,2(t*t,)2-127.8(t4T )+421)) 
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+ constant + 
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NS g(t,t, t, ) «f( o» tO t tO 04 L) (3) 


where, 


time vehicle arrives at the tunnel 


(t 
li 


=(+l. t <(t+t.) < t 
s > ^"start- 1 


-1, t < (t*ti) < tond’ 


Chapter 3 discusses ways of estimating a and 3, 
Assuming for now that this has been done, there are at least 


two ways of finding the value of ti which minimizes T. 
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l. By assuming the form of £(M t,t,)and g(N,t,t,); 
it is then possible to solve Eq. 3 numerically for the opti- 
mum value of ; An example of this is done in Chapter 4, 

2. By computer simulation; the case for this method 
is strong. (3,20,16) The tunnel-sink system is extremely 
comple: and it is impossible to describe mathematically 
exactly what will happen to a vehicle which enters the 
system, With simulation it is not necessary to guess f and 
g, nor is it necessary to assume first in, first out, Asa 
vehicle enters the tunnel, its expected time in the tunnel 
is calculated, When this predicted time is reached, the 
vehicle is simply transferred out of the tunnel into the net- 
work of roads leading to the sink and its time of arrival at 
the sink is predicted. It has been assumed that a decision 
is made for only one vehicle at a time. With simulation it 
is possible to determine the effect of letting any number of 
vehicles through to the same sink, Time of arrival at the 
tunnel (t) is really a random variable. With simulation, 
the modeler can study the effect of varlations in t as well 
as Variations in the initial conditions in the tunnel and in 
the sinks, Such a program has been written. The program 


Listing is in Appendix A and the logic flow chart is in 


Appendix B. 
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Before the experimenter can apply the theories of 

this chapter to a particular system, he must be satisfied 
that he can predict the expected waiting time in the traffic 


sinks, This problem is discussed in the next chapter, 








CHAPTER 3 
DRAWING SAMPLES FROM THE TRAFFIC SINKS 


3.1.  xponential Model 

As was discussed in Chapter 2, it is assumed that 
the sum of waiting time in a traffic sink, servicing time 
and exit time (w) is equal to ae", For any specific 
traffic sink, the experimenter must be able to estimate a 
and 8 and obtain some idea of hov well his model represente 
the conditions of the physical sink, The purpose of this 
chapter, then, is to present a sampling technique by which 
these aims can be achieved, 

To be specific, the model for total time in the sink 
is: 

m X ae "ci where € ,vN(CO, 57), O en m. 
andm = the number of observations, 
It is assumed that; 

a, deviations arenormally distributed, 

b, deviations have a common variance 62. and 

c. deviations are independent, 

The actual graph of w vs N depends, of course, upon 
the specific values of a and B, but in general,it can be 


represented as in Fig. 9. 
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Time in Traffic Sink 
(Minutes) 
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FIGURE 9. 
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Time in Traffic Sink (w) vs Number of Vehicles 


in Traffic Sink (N) 


Since N(0,6°%) e „EM exp (blues) a2) the Likee 


[IT 
| a d 
lihood functior is Ls(L)m/ 20 dm. exp(=". $ (w «ue 
2n 6 isl' i 


NIS e 
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Since *, and W, are observable values, there are three para- 


meters which must be estim^ted: a and 8, and 5^. 
mized by selecting a and 8 such that: 
m m 
BN; ^ 28N 
2 y e Puer a Se O 
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and 


L is maxi- 
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Solvins Eqs. 4 and 5 simultaneously yields the 
maximum Likelihood estimates of a and 8. It is interesting 
to note that these are identical to the least-squares 
estimates of a and 8. The sum of the squared deviations of 
the observed values of w, from the predicted values is 
m 2N 4,2 
S(a,6) = E, 1), 


tivies of S with respect of a and 8 and equating them to 


(w,-ce Upon taking the partial derivi- 
zero to minimize S, Eqs. 4 and 3 again result, 

To obtain the maximum likelihood estimate of 5^ one 
simply maximizes the partial of L with respect of 52, The 
result is Be z Ka (v, pe 1)? fm where m is the number of 
(w,,N.) two-tuples observed and ' over a parameter represents 
an estimate of that parameter, It is not possible to deter- 
mine whether or not Š is biased; this requires a knowledge 
of the bias in a and d. However, maximum Likelihood 
estimators are, in general, not unbiased, 

One advantage of maximum likelihood estimators is 
that, under fairly general conditions, it can be shown that 
they are efficient estimators, (26) From this it follows 


that maximum Likelikhood estimators are also consistent, 


3.2. Determination of Sample Size 
At this point it would be desirable to have a backe 


log of samples from which it would be possible to obtain an 


advance estimate of experimental error (82), From this the 
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experimentar could determine the sample size required to 
satisfy his precision requirements, If such a backlog of 
samples is not available; it will be necessary to obtain an 
advance estimate of 52 in another way. 

The proposed methods of drawing samples and estimate 
ing parameters will be demonstrated by means of an example, 
The traffic sink selected was a Circle E food market in 
Tucson, Arizona, It is not hypothesized that the Circle i: 
represents a New York traffic sink; however, the food mare 
ket does possess some qualitities similar to those of a 
traffic sink, rirst, most people enter Circle E with 
similar missions -- to buy one or two items, In New York, 
trucks entering traffic sinks have similar missions also. 
Second, there is reason to belleve that the congestion 
created by the presence of many cars in a Circle K parke 
ing lot will cause the rate of change of waiting time to 
increase with N, 

ror every entering vehicle, it is necessary to 
determine w, and N,, A convenience sample of size meS was 
drawn and the results are shown in Table 1, 

A modified steepest descent method using an iterae 
tive search was used to solve the nonlinear equations l and 
2. The following results were obtained: 


a = 0.726, 8 = 0,64 and 52 = 0,190, 





TABLE 1, Results of Pilot Semple 
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Now, using this estimate of experimental error, it 
is possible to determine the sample size necessary to satis- 
fy the precision requirements, These requirements may be 
worded in several ways. <A few examples will be given, 

3.2.1. Requirement I: 

Determine the sample size (m), such that E(w) can be 


estimated within % minute 99,5% of the time. 
' 


1 


Since (T 995. um d = % minute, 
(t.995 m=1)2(32) 2 
uu a? "LE 995 mer? (0-19). 


It is necessary to iterate on m until the equality 
is satisfied, It was found that with m-2, (t gos 7) (.076) 
equals 9,75, with m=9, St 995 g) 60:76) = 8,51, Therefore, 
for these precision requirements a Sample size of 9 is 


necessary. 





2 
3,2.2 Requirement II: 
Determine the sample size such that C = 99% confi- 
dence limits on E(w) will have, with probability p=,99, a 
length less than or equal to one minute, 
¡ere we use the relationship developed by Harris, 


Harvitz, and Mood; i9) 


dec Gare et Ein mel vy) 


where 
v+l = pilot sample size = 7, and 
d z l minute, 


E Y 
Fence, D = (font men? (Fon met 7) 1 
m 


Again it is necessary to iterate an m until the equality is 
satisfied, In this example, a sample size of 35 will fulfill 
the precision requirements, 

3.2.23. Requirement III: 

Determine the sample size required for the probabil- 
ity of a 90% confidence interval on E(w) having a Length 
d<8/5 to be equal to 95%, 

Letting E be the estimate of £(w,) it follows that: 


= ^ t - 
DL w4-E(w41)) < — EE = 0.90 and 


Jm 
P( w,-E(w,)) < 6/10) = 0,95, 
Rence, 


er 6/10. 
Jm 
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Io find the minimum value of m that satisfies the given 


requirements, set Szin. = 6/10, 


Squaring both sides Tn sinplying, Eq. 6 is developed: 


ns? m2 
— Be nn 
6 100t^g 95 mel’ (6) 


2 
It is now observed that ag v^ 1° 
Š me 


For values of v greater than 30, it is possible to 
use the fact that TE - /2V-T)2, is distributed 
epproximately as N(C,1) to obtain the following approximation, 
X = TT? 
Now, incorporating this result into 2 6 the result is, 


BE. = %(Z gs */2v-1)^- 


SE .95 mel 9 


x = s o5*/2m-3)* « Ri o" 
-95,mel 100t .95 „mel , 
hence 


= 50t¢ + ne 3)2. 


. 95 mel e 
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For values of m greater than 30, the approximation 


t 95,m-1'V ^.95 is valid. 


Therefore, m^ = 502% 9s(z 95 +/2m0 3) 2 and because a negative 


sample size is meaningless, 


m = M 2.95 (Z os + J/2me 3) , 
m2 - 2(19.15)m + (19,15)? = (11,62)2(2m=3), 


n^ io mn 3775 +0. 


3l 





Therefore: 
m = +310. 33/310. 32-4(1)(775) 
? ; 
a = 310.3 + 303 


2 o 
m = 305 or 4. Upon substituting these two sample sizes 
into the precision requirements it is found that m= 4 is 
too small and m = 305 is the only acceptable solution. 

It is readily seen that Requirement III is by far 
the most restrictive; it does not require an advance esti- 
mate of 52. It is felt, however, that few people would word 
their precision requirements in this manner. Therefore, for 
this example, Requirement II was selected and a systematic 


sample with a random start of size m = 35 was drawn, Ths 


results are summerized in Table 2, 


TABLI 2, Results of Sample 
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3.3. Estimation of Parameters 


eege. == ar 





It is convenitent to use a and 8 from the pilot 
sample as initial guesses in the iterative technique for 
solving Las. 4 and 5 for the final estimates of a and 8, 
The results obtained were: 

a = 0,591 minutes, B 2490.390, SS = 0 599 minutes 


9 


1 
6 = 0,774 minutes, 


A = 0,591 minutes, Wo = 1,066 minutes, 
wz = 1.922 minutes, v, = 3,468 minutes, 
We = 6,256 minutes, 


The experimenter is now in a position to test assumption c 
in section 3,1. that the deviations of w, at each value of 
N; have a common variances 62, For this the author will use 
Bartlett's test of Homogenelty of variana, |) 
mee this test, 

k/(1-L) *V /(X^k-1 where 

K 2, 3026 Oz df,)1log so - (cdf, log 3X2 

ee : =+ 

ADA. 


Kk = the number of populations of w's, 


am si , E + 885 + eee + SS, 


D zdf; 8 


In this example, 
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_ 1,078-0,033 - 0 08 
== SC DID GE, 


s = 8,97/30 = 0,299, 


los s^ = -0,324, 
p 
zdf, log Si = LO LESS 


K = 2.3026 (=15.72 + 16,15) = 0,99, and 
X? = .99/1.08 = 0.917. 


E 
2 


Since X? 9S 4" 9.48773 one cannot reject HO: M = ô 


coe = Be and the assumption seems to be valid. 
It is also possible to calculate 95% confidence 


limits on Wi. 


SACH - 367 


L 


' 
(5). 7 
T IDEE 
1 va 


ror the Circle E data, 


SH eee aro 5, 
d, = 0,48, d, * 0.56, and 
d. SOL. 


Since the experimenter is using the same estimate of 
experimental error for each array, the length of the confi- 
dence interval depends only upon ni. The graph of Vi vs Na 


is shown in Fig, 10 along with the 95% confidence at each Wi. 
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3.4. ¿Evaluation of the rodel. 

The author has estimated a, 3 and ô for the traffic 
sink, It is now desired to examine the goodness of fit of 
the model to the physical system, <=xperimental error (6) is 
always an indication of goodness of fit, The value 5 =0,774 
minutes may be pood, it is impossible to say how good how= 
ever , since there is nothing to compare this with. As was 
stated earlier, what is really hypothesized when it is 
assumed that z(w, ) = aePN 4 is that the rate of change of w 
(as well as w itself) increases as N increases, This assump- 
tion can be tested by comparing the model with a Linear 
regression model, 

The Linear model is wy 3) =a + BN; + Di +E ij 


j th th 


where Wi is the observation from the i" population, 


D, is the deviation from linearity of tho ¡ta 
population, and C, is the experimental error of the 1 
population, Again it is assumed that the experimental error 
at each value of N, is normally distributed with the common 
variance 52, This model has the advantage of testing 
directly the hypothesis that the deviations from linearity 
are equal to zero, (18) ihe analysis of variance for this 
model, using the same data is given in Table 3, 

Since E (MS departure from linearity) = + (k-2)"1+ É BIIN 


i=] ¡2 1 


and E(MS Residual) = X we can test HO:D =D,=* * *=D=0 with 


1 
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TABLE 3. Analysis of Variance for Linear Model 
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Regression m, (N,-N) 252.26 
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the ratio 


T + LS Departure from linearity 
MS Residual ° 


- = t t 
The results are N = 1,69, w = 2.56, 6 = 0,92, 6 = 0.97, 
W = 1.25, w, = 1.87, wW} = 2.33, w, = 3.79 and ws = 5.67. 


Y - : 
F = 6.395 = 5,75, and J GINE cT - 4.1/1. Hence HO is 


rejected, and the author concludes that the linear model is 
not adequate for predicting E(w,). It has not been shown 
that E(wi) = ae”Ni is the best possible model for a traffic 
sink of the type considered, But it can be concluded that 
it is better than the linear model, The linear regression 


line is shown in Fig. 10. 
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If the exponential model had not given a better 
prediction of E(w,) than the linear model, the experimenter 
would be forced to change the model of the traffic sink. 

The author feels the polynomial regression would do 
well in this situation since the experimenter has the 
ability to step the degree of polynomial until he finds 
the best fit of the polynomial models, 

At this point, it can be concluded that the expected 
waiting time of vehicles entering the traffic sink in 


question can be represented adequately by E(w, ) = der. : 





CHAPTER 4 
CONTROLLING TRAFFIC SINK CONDITIONS 


4.1. Introduction 

After estimating a and 8 and satisfying himself that 
he can predict adequately waiting times for a particular 
sink, the experimenter still must answer one more question 
before the philosophy of Chapter 2 can be applied to a 
tunnel-sink system, Let Yj represent the expected nunber 
of arrivals/minute at sink, of vehicles which have traveled 


J 


through the tunnel, This value y can be estimated from 


ale 
the origin and destination data available for the tunnel and 


must be the averaze over a certain time interval within 


start’ end 
of arrivals per minute from all sources other than the tunnel 


(t E .). Again letting R equal the expected number 


for sink ; the question is then, does the tunnel supply a 
large enough proportion Y GO 9 j), of the expected sink 
traffic to ensure that traffic controls at the tunnel vill 
be reflected in the sink? 

If 9075 of the traffic in the sink enters via the 
tunnel, then the experinenter knows he has sufficient con- 
trol over the sink to employ the method of Chapter 1, lowe 
ever, if only 1% of the sink traffic comes through the 
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tunnel the experimenter can do almost anything to the tunnel 
without affecting conditions in the sink, In other words, 
there is a critical value A* of A above which the experi- 
menter does not have sufficient control of the sink condi- 
tions to nerit controlling the tunnel traffic. For 1>1*, 
it is expected that the optimum waiting time in New Jersoy 
will always be zero, This critical value of À must be 
estimated separately for every sink and depends upon a, 8, 
the expected operating conditions of the sink, and the 


expected rush hour conditions in the tunnel, 


4.2 Estimation of the Critical Value of À 
In order to estimate A*, it will be assumed that 


Eost, t), the function which predicts the arrival of 
vehicles which have left the tunnel but have not yet arriv- 
ed at the sink, and g(ii,t,t,), the function which predicts 
the dissipation of vehicles initially in the sink, are 
known. Again let M be the nunber of vehicles which have 
enterad the tunnel but have not yet arrived at the sink and 
Let N. be the number of vehicles presently in the sink, 

A reasonable first approximation for f and g is to assume 
that they decrease linearly from their initial values, If 
N is depleted linearly and reaches zero in U^ minutes and 


Ng decreases linearly to zero in Ui minutes then: 
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1,t,<U, 


SÉ, E pa > WE eee LLL 
in 


u = aN ei? 1 0/7 
Eli , 8,41) d og 01/99) (9/1 )k, 


where 0 = explected flow for the next U, min., 
m = expected no. of vehicles in the tunnel for next 
ü. min., and 


9 


Ga 


"Zr 


Consider the second equation. The term (M erit, /U5) simply 
estimates the number of vehicles bound for the sink which 
have not yet arrived, If t, = 0, this term is p and if t,- 
U, this term is zero, Since t3 is assumed to be a constant, 
vehicles enter the traffic sink with the same flow with 
which they exited the tunnel, If all Nr vehicles were orice 
inally in the tunnel, then (it e. 64/05) /M indicates the por- 
centage of vehicles in the tunnel with this sink as their 
destination, Since @ vehicles/minute exit the tunnel, then 
(iig eris t1 /U5) (8/M) vehicles exit the tunnel each minute for 
this particular traffic sink, and ta minutes later they 
enter the sink, Then, in general, 


_ 3 Liostoltistä) , y „N t : 
w = ae pa (9/F)k^*(N, a eg), 
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ihe tern t4 is introduced to account for the time delay 
between exitins the tunnel and entering the traffic sink. 
Notice thet the vehicles which enter from sources other than 
the tunnel are never removed from the sink, this implies 
that, besides the firstein, first out assumption, it must be 
assumed that these vehicles will not depart from the sink 
before ti reaches its maximum values,  Althoush these 
assumptions seem reasonablo, they are not strictly necessary 


as will be shown later. 


6.3. ¡xample of ;stimation Procedures 





To demonstrate the method of determining A^, it is 
hypothesized that a scaled down version of the Lincoln 
Tunnel supplios the Circle K with vehicles. The Circle ii is 


used as a, B and 52 


have already been estimated for it, It 
is assumed that t, = l minute and the upper limit on t; isg 
five minutes, Figs, 11 shows the plot of E(w) vs ti for 
various values of A, Tine in the tunnel (t5) is sinply a 
scaled down version of Fis, 6, without shock waves. ¡he 
values of t, which minimize T are shown in Table 4, For 

A z 1.4, the optimal waiting time (t,) is always zero. This 
is because for \>1l.4 the experimenter cannot exert suffi- 
cient control over the traffic sink to make waiting worth- 


while, Uncontrollable sources overshadow all the 


experimenters efforts to reduce the number in the traffic 
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FIGURE 11 


Time in Traffic Sink (w) vs Waiting Time in 
New Jersey (ti) 


sink, Again using Fig. 6 as a guide, a similar analysis 
was made assuming that a shock wave occurred in the tunnel, 
The results were the samo and indicated A* = 1,4, 

In order to circumvent the first-in-first-out assump- 
tion and the assumption that no vehicle which enters froma 
source other than the tunnel will exit the sink before 
the method of depleting N must be modified, It 


t = 
1 "imax? 
can be assumed that a constant percentage of the vehicles 
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TABLE 4, Linear Depletion of f and g 





ti optimum T minimum (minutes.J| 
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which are in tho sink at the beginning of one minute will 
have erited by the end of that minute, At the beginning of 
each minute, vehicles enter the traffic sink from uncon- 
trollable sources and possibly from the tunnel and combine 
with vehicles already in the sink, By the end of this 
minute a certain percent of the vehicles have departed. 
Besides the previously mentioned advantages, this approach 
does not stipulate the manner in which vehicles leave dur- 
ing a riven minute, This method does require an updating of 
N at the end of each minute before vehicles enter for the 
next minute, but still falls under the category of assuming 


the form of f and g. 
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Tig 12 shows w as a function of Ei if no shock waves 


are encountered in the tunnel and if the constant percent 
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Time in Traffic Sink (w) vs Waiting Time in 
New Jersey (tj) 

depletion method is used to update N. For dA 20,6 it is 
seen (Table 5) that the experimenter does not have suffi- 
cient control over the sink traffic to make waiting before 
entering the tunnel profitable. Similar results were 
observed with a shock wave in the tunnel, T minimum was 
increased but \* remained 0,6, 

;ne experimenter must now decide, from his know- 


ledge of the traffic sink, which method to employ in 
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TABLE 5, Constant Percentage Depletion 
of f and g 





t, optimum T minimum (minutes) 
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estimating A*, Tor example, the author feels that the first 
in, first out assumption is not unreasonable for the Circle 
& sampled, However, he does not believe that the first 
method can be realistically applied. While making a decie 
sion for one vehicle, some of the contontrollable vehicles 
will enter and exit while 0 < t; S El mag, Although this 
will not be true if w is expected to be larger than vico D 
the author feels that for the Circle K, A* = 0.6 is a more 
reasonable limit than A* = 1,4, 

It is not difficult to estimate A. While sampling to 
estimate a and B, the easiest method is to record the 
initial value of N and then mark the time of each arrival 
and departure, for this X can be determined, Precision 


requirements on A may necessitate drawing a larger sample 





than is indicated in Chapter 2, however no theoretical 
problems are involved, 

If X < A*, the experimenter can feel that he has 
enough control over the sink traffic to be able to reduce 


T by controlling Ei, 








CHAPTER 5 
SUSMPARY AND CONCLUSIONS 


29.1. Summary 

Ihe hypothesized model for predicting total time in 
the system is T - €, + t. + constant + ae’ N(t1), As Ty 
increases, N and hence E(w) decrease, The objective is 
to find the value of ti for which the reduction in T 
attributal to descreases in E(w) is offset by the increases 


in T caused by the zrowth of t, and t,(t)). 


1 

Before attemptins to use this model on an actual 
tunnelesink system, two basic questions must be answered, 

l. Can the total time in the traffic sink be 
predicted adequately? To answer this question, the experi. 
menter must sample the sinks to estimate the parameters a 
and B, using, for example, proceedures outlined in Chapter 3, 

2, Does the tunnel, over which the experimenter has 
control, supply a large enough percentage of the input to 
the sink so that the reduction of N will overcome the 
increases in T caused by ti and to? This problem is dis- 
cussed in Chapter 4 and involves estimating A and A*, 

For the tunnel, it was found that t5 increases with 
ti as in Fig. 6 of Chapter 2. For the Circle K it was found 
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t 
ta = 0.59, p = 0.591, = 0.67, and X* = 0,6, Thus, if 


a scaled dom version of the tunnel supplied vehlcles to 


the Clrcle K, the optimum t, would be 4 minutes and would 


result in an expected total time in the system of 13,3 


minutes, 
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These results are summarized in Fig, 13, 
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FIGURE 13. 


System (T) vs Waiting Time in 


New Jersey (t,) 


Since it is reasonable to believe that t4 will actually 


decrease as t 


1 


increases and thus help decrease T, the 


assumption t, = constant is a conservative one, 
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2.2.  Arplication of the odel to the Lincoln iunnel System 
ihere are several stumbling blocks which must be 
overcome bafore the idea can be applied successfully to the 
Lincoln Tunnel e New York City system, 

l. There are no vehicle waiting areas aVailable in 
New Jersey, 

2. As of this writing, the authorities had not yet 
set up a traffic control in the tunnel so as to avoid shock 
waves, 

3. Because of the sizes of the traffic sinks, the 
number of vehícles involved, and the large number of sink 
entrances and exits, the sampling problems in New York are 
expected to be very difficult, much more so than for the 
example developed in Chapter 3, 

Due to the myriad of entrances into any New York 
sink the question of controllability posed in Chapter 4 is 
very real. It is unlikely that the Lincoln Tunnel supplies 
a large enough percentage of the traffic into the Central 
Business district of New York City to make this method of 


control applicable to all 24 traffic sinks, 


5.3. Extension to U, S, Naval Ships 
The theory and methods discussed in Chapter 2 apply 


to any transportation system in which controlling the 


vehicles as they enter will affect the conditions within 





50 
the system that these vehicles will encounter, ‘to the 
author, one of the most obvious extensions of these control 
procedures is to Naval ships entering port through the cone 
straints of a channel. Many of the problems encountered in 
New York are alleviated in this channeleharbor system and 
some new difficulties arise, but the basic control theories 
remain the same, 

The traffic sinks within the harbor are the deep 
draft piers, the shallow draft piers, fueling piers, 
anchorages, possibly ammunition loading areas and possibly 
several points within a shipyard. There are fewer sinks 
than in New York and they are more distinct and easier to 
define, The sampling procedures are simpler than those in 
New York and probably would require Little more than sample 
ing the watch officers’ logs of selected ships, 

There would be no requirement to construct waiting 
areas at the channel entrance, and the problem of controle 
lability is completely removed. The Senior Officer Present 
Afloat (SOPA) has absolute control over all vehicles which 
would enter the traffic sinks, Thus A would be zero and A* 
noneexistent. Because of the speeds and spacings involved 
it is much easier to predict conditions in a channel than 
to predict conditions in a tunnel, 

The interactions between ships attempting to be 


serviced (to moor) in the same sink are just as real and as 
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bothersome as the delays caused by trucks trying to off load 
cargo in the same area, It is believed that because of the 
simplifications inherent in a channel-harbor system, the 
control procedures outlined in the preceding chapters can 


be applied successfully to these systems at this time. 


5.6, Conclusions 

The model of Chapter 2 is by no means complete. Tor 
example, no mention has been made of the Holland Tunnel, 
another facility under the control of the Port of New York 
Authority, 

The initial study of a large problem, must by 
necessity, involve many simplifying assumptions. These 
enable the experimenter to isolate major problem areas as 
was done in Chapters 3 and 4, 

In conclusion, it is believed that the idea of 
reducing traffic delays in New York City by controlling 
vehicles before they enter the tunnel is possible and 
worthwhile, Authority officials feel that it will be years 
before they actually begin experimenting with the ideas 
developed in this thesis, But the problem of traffic cone 
gestion , outlined in Chapters 1 and 2, must someday be 
faced, It is hoped that this study constitutes a first 


step, 
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